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ABSTRACT 


The design of airplanes, ships, automobiles, and so-called "sculptured 
parts" involves the design, delineation, and mathematical description of 
bounding surfaces. A method is described which makes possible the description 
of free-form doubly curved surfaces of a very general kind. An extension of 


these ideas to hyper-surfaces in higher dimensional spaces is also indicated. 


This surface technique has been specifically devised for use in the 
Computer-Aided Design Project at M.I.T., and has already been successfully 


implemented here and elsewhere, 
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SECTION I 


INTRODUCTION 


The purpose of this work is to present the mathematics of a certain class 
of surfaces which are suitable for the design and description of arbitrary shapes. 
In the past, the subject of surface mathematics has been investigated, in ana- 
lytical geometry and in differential geometry, from the standpoint of the analysis 
of geometric properties of surfaces that already exist, but very little literature 
has been produced on the subject of the creation of such surfaces. As a typical 
example, the design of the hull of a racing yacht requires the description of a 
surface of considerable subtlety and complexity, and the process is traditionally 
carried out by purely graphical procedures which are exceedingly laborious, 
since they entail a large amount of trial and error iteration in order to assure 
that the surface is completely described, and is smooth and "fair."' The design 
of automobile bodies and airplane fuselages is similarly tedious and time con- 
suming, although mathematical techniques have been applied to aircraft design 


for a number of years. 


A few papers have been written on the subject of fitting existing ship hull 
shapes by means of various types of polynominals, with the two-fold purpose 
of smoothing and interpolating the information contained in preliminary graphi- 
cally derived hull lines, and of replacing this graphical information with for- 
mulas and equations that will permit further analytical techniques to be applied, 
such as structural analysis and the discipline of fluid mechanics. But these 
mathematical techniques are applicable only when the surface has already been 
designed to some degree of completeness, so as to furnish enough information 


to make the mathematics work, 


The mathematical structure of the surfaces to be described in the follow- 
ing discussion has been devised to implement the surface design process itself, 
so as to make it, from the designer's standpoint, extremely natural and easy. 
The designer himself need not know or care about these internal mathematical 
details, any more than he needs to know the specific composition of the pencils 


with which he writes or the mechanics of the splines with which he now draws 
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curves. The mathematics is relatively simple, but it is nevertheless too com- 


plicated for hand calculation, and is designed for use on a computer. 


In the design of a three-dimensional object , whether it be an airplane fuse- 
lage, an automobile body, a ship's hull, or a single sculptured part of a machine, 
the designer requires a system which will permit him to define a surface with a 
minimum of input information, and then to modify this surface, if he feels so 
inclined, either by changing the original input, or by adding more design con- 
straints to the system. 


As a specific example, suppose a designer wishes to design an airplane 
fuselage, using the SKETCHPAD system. A228 He would like to be able to draw 
the outline of the airplane as seen from the side, the outline of the airplane as 
seen from above, and some arbitrarily selected section midships. With these 
three arbitrary curves designed, he would like to have the computer automati- 
cally and immediately generate a "fair" surface and display this surface to him 
in sufficient detail so that he could make appropriate judgments. If the surface 
so generated does not satisfy him, he would perhaps like to modify his original 
design curves, or else he might perhaps like to add other new sections and have 
the computer automatically and instantly re-fair the surface to fit this additional 


information. 


The following sections describe a very simple, flexible and general class 
of surfaces which are able to fulfill these requirements. It will be shown that 
a single algorithmic structure and essentially only two symbol types serve to 
provide the following features: 


1. Smooth, fair surfaces can be defined by a minimum number of curves, 
and then adjacent surfaces can be designed to match position, slope, 
curvature, and indeed any desired order of derivative along the ad- 
joining boundaries. 


2. The design curves that define the surface can be of any kind whatso- 
ever, including circles, second-degree curves, polynominals, trans- 
cendentals, and also sketched curves with no known mathematical 
formula whatsoever. 


3. Some classic surfaces are not necessarily members of the family 
of surfaces to be described; nevertheless, these classic surfaces 
can be matched along their boundaries to any order of derivative 
desired. 
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4. The arithmetic involved in constructing these surfaces is extremely 
simple and, we have found, easy to implement on a digital computer. 
It also lends itself to special-purpose computing hardware, such as 
digital or analog differential analysers. In addition, by virtue of the 
form of the algorithm, the parameters that define the shapes are ex- 
tremely easy to compute. (In some cases they may require no com- 
putation at all.) 

We intend to develop a method to construct complex arbitrary surfaces by 
piecing tegether surface "patches.'' Each such patch will be defined by four 
boundary curves, in principle, although it is harmless for one of the boundary 
curves to be degenerate, and to appear as a point instead of a curve segment. 
In the design of a surface, it is intended that the designer begin with a single 
surface patch, or a very small number of patches, and then subdivide these 
regions with additional design curves defining boundaries of smaller patches 
only when the internal surface needs modification. This is somewhat at variance 
with the customary procedure for mathematical curve fitting and surface fitting 
of existing curves and surfaces, in which a relatively large number of surface 
points already defined by some other procedure are used to obtain mathematical 
expressions for a surface that best fits them. Instead, the system to be de- 
scribed is intended to be used by the designer at the outset, in the process of 
designing the surface, rather than later on as a means for making it mathe- 


matical. 


This is not to say that the surface-patch technique cannot be used to for- 
mulate patch-wise mathematical expressions for existing surfaces, but rather 
to indicate that the primary purpose of this surface technique is to facilitate 
the initial design process itself. 


When the design process is completed, the surface will be completely 
mathematically defined, since this definition occurs automatically and concur- 


rently with design. 


Ordinarily a ship's hull or an airplane fuselage is described by certain 
important curves such as, in the case of the hull, a keel curve, a midships 
section, and a curve representing the sheer or deck line; these curves are 
sufficient to determine a surface, since they form the boundaries of a surface 
patch. However, ordinarily this primary surface will not have certain desired 
characteristics, and it will have to be modified by introducing additional 
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information, such as for instance one or two other section curves. When these 
additional curves are introduced, the surface algorithm permits the computer. 


to "re-fair" the original hull form to contain these curves. 


Similarly, an airplane fuselage can be designed by drawing a profile curve, 
a maximum half-breadth curve, and a mid-section of the fuselage. Again these 
curves suffice to define a primary surface, which in most cases will require 
modification by the addition of a few more curves to make more explicit the 
designer's wishes. As these additional curves are introduced, the original 
surface will be sub-divided into patches, but the algorithm will automatically 
insure continuity of surface slope and curvature (if desired) and will incorporate 
these additional curves into the surface automatically. This should make initial 
surface design virtually painless, and is intended to remove the tedious process 
of surface fairing as it is now practiced in naval architecture. Airplane fuse- 
lages are usually somewhat simpler shapes than, say, yacht hulls, and fora 
number of years second-degree curves have been successfully used for fuselage 
design; on the other hand, naval architects have steadily resisted the use of such 
methods in their work, since the complexity of yacht shapes makes it necessary 
to pay attention to the irksome details of the geometry involved, and second- 


degree curves prove to be cumbersome in such applications. 


The system that is described in this report is’ intended to furnish the flexi- 
bility that second-degree curve techniques lack, and to remove almost entirely 
the need for the designer to be an analytical geometer. With this system imple- 
mented on a computer, there is reason to believe that the computer can take 
over all of the geometrical and mathematical burden of the design process, and 
leave the user free to be a sculptor assisted by an exquisitely skillful mechan- 


ical slave. 


Ultimately, when a graphical input-output hardware for a computer is 
available in the engineering design office, these methods will permit designers 
to delineate complex shapes with great ease, by simply drawing the salient 
curves that define and describe them. Already experiments along these lines 
are in progress in a few isolated laboratories both in universities and in 
industry. Very soon the two severe handicaps that have inhibited the wider 
use of such graphical devices will be removed. These inhibiting factors have 
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been high cost for the terminal hardware and small size of the working area. 
Rapid strides are being made on both these fronts, and within a few years it will 
be possible not only to draw on a virtually unlimited drawing surface, but to 
draw objects directly in three-dimensional space, and to view these constructed 


objects as one would view an actual physical thing. 
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NOTATION 


We shall in what follows relate the x, y, and z coordinates of points on a 
surface to two independent variables u and w, so that we could write 


x =f (u, w) 
y = gu, w) 
z = hu, w). 


If the functions f, g and h were specified, then for a pair of values of u and w, 
a point in space would be defined. If we held one of the independent variables 
fixed, say w, then by allowing u to vary, the point in space would trace out a 
curve. If subsequently we set w to a new fixed value and again allowed u to 
vary, we would trace out another curve, andsoon. Clearly by stepping the 
values of w by small increments and allowing u to vary after each such step, 
we could produce a family of space curves that would lie on the surface and 
define it. All that is needed is some convenient and systematic way of arriving 
at the functions f, g, and h. 


It will turn out that the form of all of these three functions is the same; 
only certain internal numerical values are different. In vector notation we can 
write 


[x y 2] = E (u,w) g (u,w) h (u,w)} 


~ --Since- V = [x-y- 2} -is-a suitable conventional abbreviation for the vector 
on the left, we introduce a similar abbreviation for the right hand side: 


ww) = [Faw gt.) b (uw) 


Here, in the abbreviated symbol on the left, we shall omit the comma between 
the two letters. Later on, when no ambiguity can arise, we shall omit the 
parentheses as well, and write simply uw to stand for the vector. It is to be 
remembered that uw does not stand for the ordinary product of the two quan- 
tities, but is merely a bi-literal symbol standing for a vector whose components 


are functions of the two variables. 
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We plan to build up surfaces by adjoining surface "patches," in an analogy 
of the piecewise fitting of complicated curves by curve segments suitably joined 
together. Accordingly, we shall at the beginning focus our attention on one such 
surface patch. To simplify arithmetic, we shall stipulate that the independent 
variables, or parameters, u and w can take on only values between 0 and 1. 
Then a surface patch can be considered to be a surface segment bounded by 


four space curves, (0 w), (1 w), (u 0) and (u 1). 


11 
01 ul 


lw 


10 
"00 ue 


Here, typically, the symbol (Ow) stands for. the vector describing the x, y, and 
z coordinates of points along the curve generated by allowing w to vary conti- 
nuously from 0 to 1, while u is held fixed and equal to 0. 


We shall introduce two scalar functions, F, and Fi each a function of a 


0 
single variable. These will be referred to as "blending functions" for reasons 


that will become clear. 


In order to compress the surface equation, and the proofs that we wish to 
demonstrate, we shall use a kind of indicial notation; we introduce the indices i 
and j, which can assume only the values 0 and 1, and we invoke the customary 
summation convention for terms with repeated indices. This convention in our 
case simply means that when an index is repeated in a term, we write out all 
the possible terms that the actual indicial values generate, and then add them. 


SECTION II 
THE SURFACE EQUATION 


With these conventions and notational peculiarities in mind, we write 


(uw) = Gw)F, (a) + (uF (w) - Gj) F, (uF i) 


j 
(Typically, the first term on the right expands as follows: 
(iw) F, (a) = (Ow)F,(u) + (1W)F, (u). 


Thus the complete expansion would consist of eight terms, if carried out.) We 
shall proceed to demonstrate that this surface equation represents a surface 


that contains the four boundary curves, and is thus defined by them. 


We must make a stipulation, a weak one, on the nature of the blending 


functions Fy and Fy: 
F (0) =1 Fy) =0 
Fi(y= F,(0) =0 


A further stipulation is that F 0 and FY be continuous and monotonic over the 


interval. 


Now set u = a, where a can only be either 0 or 1. Then, substituting in 


the surface equation, 
(aw) = (iw) F(a) Ba (aj) F(w) * CE OYE 9) « 
Consider F, (a) which occurs twice in the equation, By the stipulation, ifi=a, 
F(a) =1, 
Otherwise, if i #a, F(a) =0. 


Hence all terms in the expansion that contain i # a vanish; we can set i = a and 


what remains is 
(aw) = (aw) F, (a) + (aiyF j~ = aj) F@)F Wm) 
= (aw) + (aj) F (w) 7 (aj)F(w) 


= (aw). 
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This shows that for a = 0 or 1, and hence (aw) = (Ow) or (lw), the surface equa- 
tion reduces to an identity. This implies that the surface contains its boundaries. 
An entirely parallel argument would show that the equation also reduces to an 
identity for the other two boundaries (u0) and (u1). 


Provided a pair of functions F 0 and Fy, are chosen once and for all that 
satisfy the stipulations, the surface equation may be constructed immediately 
and uniquely for any set of boundary curves (u0) (ul) (Ow) and (lw). It is to be 
observed that no restrictions have been placed on the form of the boundary 
curves; there is perhaps the restriction that they form a closed boundary, at 
least at the corners (ij) = (00), (01), (10), and (11) otherwise there will be mul- 
tiple values within the surface segment; similarly they should be continuous 
functions, but apart from these rather obvious restrictions, they can be of any 
shape whatever, including curves that can only be represented by tables of 


values. 


We can gain intuitive insight into the nature of such a surface if we look 
at one of the terms, say uF, (Ww). 


We have the expansion 
(uF (w) = (U0)Fo(w) + (ul) F, (w). 


This represents a weighted average of the quantities (u0) and (ul). When 
w=0, F 0) = 1 and F (0) * 0, and the expression becomes simply (u0). As w 
increases, the weight of Fy) decreases, while that of F, (w) increases, so that 
the surface partakes of the nature of both boundary curves. As w approaches 
the value 1, the influence of (u0) on the shape of the surface gradually disappears, 
while the influence of (ul) gradually becomes dominant. Finally, at w =1, the 
curve (ul) represents the shape of the surface. We can say that the surface is 
generated by a gradual transition from (u0) to (ul), and that these two curve 
shapes are "blended" together by virtue of the blending functions F 0 and Fi 
This discussion is somewhat oversimplified, since we have omitted the term 
(iw) F, (u) and it too plays a part in determining the shape of the internal surface, 
as does of course the term involving the corner coordinates, (ij)F; mF (w). 
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The entire surface equation is seen to be symmetric in u and w, and by 
virtue of this and a secondary symmetry in the functions Fo and F, , we can 
abbreviate proofs about the behavior of the surface along all boundaries by 
exhibiting a typical proof for any one boundary. 


3.1 BOUNDARY SLOPE CONTINUITY 


It is our aim to design and delineate complicated surfaces by adjoining 
surface patches, in a plecewise fashion, Consider two such patches A and B, 


with a common boundary. For patch A the boundary is (1w); for patch B it is 


(Ow), and the vectors of coordinates are equal, 
A (lw) = B (Ow). 


Then the two patches will be continuous across their common boundary. They 
will however in general be discontinuous in slope across the boundary, and we 
wish to investigate this and make some amendments that will correct this dis- 


continuity of slope, 


We take the partial derivative with respect to u: Our symbolism for this 
partial derivative is (uw), = a) , and when we substitute, say, u=0, we can 
write (Ow), to mean the value of the partial derivative so obtained, Then 


(uw), = Gw) Fi (u) + (ad, F(w) ~ i) F} (WF (w). 


j 


Now substitute u =a =0 or 1, as before. 


(aw),, = Aw) Fi (a) + (ai), Fim) - (1) Fa) F 0), 
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If we now place additional constraints on the blending functions, that their first 


derivatives 
F(a) =0 (a = either 0 or 1) 
we obtain the result 


(w), 


(ew), = (2)),F 
all other terms vanishing, 
This implies, for example, that when a = 0, 
(Ow), = (00) Fo(w) + (OF (w), 
or, the derivative anywhere along the boundary in the u direction (across the 
boundary) depends only upon the derivatives at the end-points of the boundary; 
it is entirely independent of the shapes of the four boundary curves, including 


¢ 


the boundary (0w) itself. a 
Thus for the two patches A and B, if 
A(10), = B(00) 
and 
A(11) = B(01), 


i.e., if the boundary curves are continuous in slope in the u direction at the ends 
of the contiguous boundary between patches, we are guaranteed to have 


A(1w) = B(Ow), everywhere along the boundary regardless of the 
shapes of the boundary curves of A and B. This is a remarkably powerful and 
useful property, achieved at the slight expense of extending the stipulations on 
the Fi. 


Similarly, the second derivative with respect to u is 
= Tt - i tf 
OW) co OY Fi (u) + MDF (ij) Fi SUE) 
and if we further stipulate that Fi (a) = 0 we obtain 


OH) = Day F(M- 
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This establishes second derivative (or curvature) continuity as an auto- 
matic and inherent property of adjacent patches, provided their boundary curves 
have this kind of continuity at the end-points of the boundary, It is easy to see 
that we may escalate in this way to any level of derivative continuity we wish 
along contiguous boundaries. 


3.2 SLOPE CORRECTION SURFACE 


The surface equation already described is very general, in the sense that 
it can contain virtually any boundary curve we wish, and it has certain benign 
properties of derivative matching along boundaries; nevertheless it is not a 
universal formula for all surfaces, and there are many that do not belong to its 
family. We have already seen that surfaces generated by the surface equation 
have a definite intrinsic slope along boundaries, whose variation is rigidly 
prescribed by a simple formula. Obviously surfaces exist whose boundary 
slopes do not match this intrinsic slope, except at the end-points of boundaries. 
Nevertheless, we wish to be able to patch together such other surfaces with our 
special surfaces, so as to have slope continuity (or continuity of any level of 


derivative). 


To do so, we introduce a new surface equation, describing a slope-correc- 
tion surface, which when added to the first surface equation has the property of 
leaving the boundaries unchanged, but causing the derivatives across boundaries 
to vary in any arbitrary way we wish, as we move along the boundary. 


The equation resembles the first form very strongly. It is 
(uw) = (iw), G,(u) + (uj), G,() = BD wFG,(w). 


Here, typically, (iw) is a function of w only, and describes the arbitrary 
variation of the derivative with respect to u as w varies, along the curve (iw), 
and similarly for the other boundaries. The vector Gi ow represents the cross 


derivatives of the four corners. Typically, 


2 
= 2 (uw) rs 
(00) w dudw bs a : 
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The functions G and G, are again blending functions or weighting func- 
tions, but they have properties different from the functions Fy and Fy. We 
stipulate 


G(0) = G, (0) = Go(2) = G, (1) = 0 
t = ' = 
G)() 1 G}Q) 1 
’ = ' = 
Gy) = 0 Gi (0) = 0 
or in the indicial notation used earlier, 


G,(a) = 0, a andi =0or 1. 
Gi (a) =0, afi. 
Gita) el, a=i. 


We need to ensure that the vectors describing the boundaries vanish 
identically, and that the vectors describing the slope variation along boundaries 
are indeed given by the equation. The proof proceeds along precisely the same 


lines we used before. First, substitute u =a. The equation becomes — 


(aw) = dw), G,(a) + (ai), G.(w) - A), G,(aG(w). 


j 
= (2) G0) 


Consider (aj), We wish to have the correction surface leave the original 
boundary vectors unchanged, and hence the boundary vectors of the correction 


surface must vanish; ie, 
(iw) = 0 
(uj) = 0 

Then the derivatives of these boundaries must also vanish; in particular, 
dw), = 0 and then Gi), =0, when w=j. Hence (aj), =0, 


Thus (aw) = 0 indicates the desired behavior of the correction surface along a 


boundary. Similarly (ua) = 0. 
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To examine the slope variation along a boundary, differentiate the 
equation with respect to u: 


(uw),, = dw), Gi (a) + (ap, Gm) 


- 4) Gy (WG 


™ . 


Now substitute u =a 
(aw), = (tw) Gi (a) + OD yj - G)) wi 9G) 


As before, G; (a) =1if and only ifa =i, 


80 we get 


(aw), = (aw) Go (@) + (Ad, G,(™) - 8), ,G) @IG,(w) 


= (BW), + BD GU") ~ BD gy GC) 


= (aw) 

This demonstrates that the surface has the slope variation along the 
boundary as required. To make use of this slope correction surface, we must 
first determine what the intrinsic slope of the surface to be corrected is, and 
then we must subtract this slope from the desired boundary slope, to yield the 
correction slopes that enter into the equation, Thus if (Ow) is the desired 
slope, and 


I (Ow), is the intrinsic slope, then 
c (Ow), will be the correction slope, 
Cc (Ow), = (Ow), -I (Ow). 
The correction slopes C (iw) and C (uj), are the four functions that enter 


into the slope correction surface, ._The desired surface is obtained by adding 


the correction surface to the first surface: 
(uw) =I (uw) + C (uw) 


where we use the symbol I (uw) to represent the surface whose boundary slope 


is being modified, 
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3.3  HIGHER-ORDER CORRECTION SURFACES 


Analogous forms may be obtained for correction of higher derivatives along 


boundaries. For second derivative correction, the surface equation is 


(uw) = (iw)  He(u) + (uj) How) - (ij) (4)H,(w). 
uu i ww j 3 


H. 
uuww i 
In this equation, the blending functions H, have the stipulations that, for a = 0 


or 1 as before, 
H,(a) = 0 
Hi (a) =0 
tt 
H, (a) =O,if#a 
Wt 
Hy (ij) =1,i=a, 


With these constraints on the Hs it is easy to arrange matters so that this 
second-order correction surface is zero everywhere on the boundary, has zero 
slopes across boundaries, and has second derivatives across boundaries speci- 
fied by GW) a and ay whatever these functions may be. The addition of this 
surface vector to a given surface vector will then provide a means for boundary 
second-derivative correction without disturbing either the boundary shapes or 
boundary slopes. 


Although we have already carried out a similar proof for slope correction, 
it might be well to exhibit once again the course of the argument. 


First, to show that the boundary vectors are zero, substitute u =a: 
(@ =Oor 1.) 


Caw) = OW), HO) + 2D FO) ~ Cyt OH Om) 
= (85) Hy (w)- 

The term 3) refers to the second derivative in the w sense at each of the 

four corners, such as, typically, (00). As in the case of slope correction, we 


must have 


(iw) = 0 along boundaries. 
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Then (iW) oy =0, GD wy = 0, and in particular (2) wy =0, so that the 
equation satisfies the boundary condition. 
For boundary slope vectors, differentiate with respect to u: 
1 ' 
(uw), = (iw) AL (u) + (8D) ratty) = @) ww (u)H,(w) 
Set u =a; 
(0), = 8D watt"): 
We wish to have the slope vectors vanish along boundaries, so typically 
(iw), = 0 for all w. 


But then GW) aw = 0 and (2) ww = 0 by taking derivatives. The order of 


differentiation is immaterial, so 
(iw) = (iw) , and finally we can conclude that 
uww wwu 


(aj) = 0; again the right and left hand sides of the equation are 


in agreement, 
Finally, we differentiate again with respect to u: 
tt tt 

Uw) = (™) Ay (a) + UH) wnat) = GD wou (OVE ts 

Set u =a; only terms in which a =i remain: 
tt 
(aw) vu = (aw) A, (@) 7 2D) wnat) = 2) wonta (a)H (=) 
= (QW) uy 


Again we have demonstrated an identity, The escalation to any level of 
boundary derivative correction vector is obvious. 


3.4 MATRIX FORM 
The surface equation 


(uw) = EW), (a) + (HF (w) ~ GF, U)F CW) 
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may be expanded directly into matrices, to yield: 


(uw) = [uo wil Fow + Fu Ful Ow 


F 1” lw 
2 [Fu Fu] 00 01] [Few 
10 «6t} |Fiw 


In this we have treated the indicial form term by term in a straight- 
forward way. We shall in what follows omit parentheses, since no misunder- 
standing can arise. Thus typically Fou is written in place of F o™) as a matter 
of convenience and economy. Similarly, typically 00 is written instead of (00); 
the reader should be reminded that this is merely a compact way of exhibiting 
the x, y, z coordinates at point (00). 


It means; 
00 = [xcoo), y(00), 2(00)| when written out completely. 
The three vector (matrix) products are equivalent to the following three products: 


(a Fuu Fu 0 ud ul 1 


0 
0 0 0 Fyw 
0 0 0 Fw 
+ {a Fou Fu] 0 0 0 1 
Ow 0 0 Fw 
lw 0 0 Fw 
+ (i Fou Fal 0 0 0 1 
0 -00 -01 ad 
0 -10 -11 Fw 


and in this form we can perform the addition, obtaining 


(uw) = F Fou Fu] 0 ud ul 1 


0 
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It is slightly more convenient to rewrite this in the equivalent form 


(uw) = - [-1 Fou Fu] 0 ud ul -1 
Ow 00 01 Fow 
lw 10 11 Fw 


so as to avoid the awkward minus signs in the 3 x 3 matrix. 


Two facts should be noted, The leading row vector in front of the matrix 
and the trailing column vector following the matrix are transposes of one 
another, but with different arguments; the matrix represents the boundary con- 
ditions of a patch, The partition Pe a is redundant, since its elements 

10.) «=i 


must agree with uj and iw for u and w equal to 0 or 1. 


We have already suggested that we can maintain slope continuity across 
boundaries by suitable stipulations on F,, and we have also already suggested 
that when desired we can adjust slopes across boundaries by a second additive 
vector with suitable stipulations on its G. We shall now investigate the com- 
bined form of the surface equation. To do so we shall prefix a symbol to the 
vector uw to indicate whether we are talking about the first surface equation, 
or the correction surface equation, and we shall omit the prefix symbol when 


we are talking about the combined form. Thus 


uw = suw + cuw, with 
suw =the primary surface 
cuw = the correction surface 


uw =the combination, 


Accordingly, using this notational convention, we will take derivatives, with 
respect to u, of the surface equation suw in order to determine its slope vector 


in the u direction, 
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1 t 1 1 
suw, =- [o Fyu Fi u| 0 sud su - 
sOw s00 s01 Fw 
slw s10 sll] Fw 
- - -1 
[-1 Fou Fu] 0 sud sul 
0 0 0 Fy¥ 
0 0 0 Fiw 


We substitute u = 0, and obtain 


sow, = [s00, s01,| 


ad 
Fiw 

Now consider, for example, 800, and the desired 00. - The symbol 800, 
refers to the slope vector at a corner; we have already seen that at corners the 
correction surface 00, =0, and so 00, = 500. This is bourne out intuitively 
by the reflection that at (ij) corners, the two crossing boundary curves com- 


pletely define the slopes there; since this is so, no correction of slope need or 


can be applied. 


Hence we should write 
sow, = [oo, 1, | Eel 
Fw 
By analogy and symmetry we can write the remaining three statements; 


slw. = fo, li, ] FW 
Fw 


sud = Fou Fu] i 
w 
sul, = [Fu Fu] ae | 
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In order to obtain a desired slope vector along any of the boundaries, we add 


the correction surface, whose equation is 


cuw = - [-1 Gyu Gu] 0 cu, cul, -1 
cow | c00 mt Gow 
clw c10 ell G,w 

u uw uw 1 


As we have already remarked, the desired surface uw is the sum of the vectors 


suw and cuw, Hence the correction slope vector, such as, typically, cud, is 


cud. = uo, - [Fou Fu] i 
Ww. 


This is an entry in the correction surface matrix. 


Now we introduce a new fact: the corner cross derivatives of the primary 
surface equation are all zero. To show this, differentiate the indicial expres~ 
sion first with respect tou, then with respect to w, and finally set u=a, w=b, 


where a and b are as usual either 0 or 1. We have 


uw, = GW)F (u) + (ui) Fim) - DF! WF (w), and 
UW gy TAM FY (Uy + (OH, FS (w) - GFL (WF! Cw). 


Evidently this expression vanishes for (uw) = (ab). This shows that the corner 
"twists", or cross derivatives, of the original surface all vanish; it is a 
peculiarity of the first fundamental surface equation. 


Hence we can assert that 


clj =ij ; This says that the desired twists at corners are 
uw 


uw 
identical with the correction surface twists, since the fundaments surface has 
no twist, We shall use this result to replace the partition c00 col 
c10 ell 
uw uw 

with, 

00 01 

uw uw 
10 11 


uw uw 
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We can rewrite the expression for cuo_ as follows: 
w 


T et 
cuo = ft -Fou -F ju] ud 


and 


cul = [1 ~F ui -F u] ul 
Ww 0 


and of course, perhaps trivially, 


0= (1 Fou Fu] | 0 


Each of these matric products represents an clement of the top row of the 


correction surface matrix. 


Since the row matrix [L arg -F ju] is common to these three products, 


it can be factored out and introduced into the matrix [-1 Gp¥ Gul to yield 


> -Fu -Fu 4 
[ [ Bae Pea Goe Ge] 
which is the same as the vector 
[-1 Fou Fyu Gu Gul 
We replace the elements of the top row of the correction surface matrix 


by the three matrices 0 i 2, and ul 
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This causes it to become a 5 x 3 matrix, and we now have the intermediate 


result, 
= ee 
cuw = - [-1 Fou Fu Guu Gul 0 wo em UB 
| 
0 00 ol 'G w| 
Ww Ww | 0 ; 
0 10 11 lGwi 
Ww Ww : 1 | 
cOW 00 Ol 
u uw Uw 
clw 10 11 
u Uw uw 


By similar procedures, we can write for the elements of the first column 


of the correction surface matrix, 


cow = [ow 00 01, | fi | 
u u u u 
~Fow 
-F_W 
CF | 
clw, = [lw wo itt Pa | 
u u Uu ud 
| 
=r 


and again trivially, perhaps, 
0 = [000] 1 
ne 
“Fw 
When we factor out the common column matrix as before, and replace each 


entry of the column matrix o | by the above expressions, 
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we obtain the complete matric expression for the correction surface: 


cuw = - [-1 Fou Fyu Guu Gu 


suw =- [-1 Fou Fyu Go Gu] 


In this bordering process, the value of the matric product is unchanged. 


Since the pre- and post-multiplicative matrices in this equation are the 
same as those of the correction surface equation, we can add the two 5 x 5 


matrices and pre~ and post-multiply by the two vectors. We shall perform, in 
fact, 


uw = suw + cuw, and obtain 


— a - -1 
uw [ 1 Fou Fu Gu Ge | 
Fw 
0 
Eye. 
Gow 
G.w 


This is a general expression for a slope~matching, slope continuous surface patch 


with entirely arbitrary boundaries and entirely arbitrary slopes across these 
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boundaries, There are no stipulations whatever on the nature of the boundary 
slope function. The stipulation on the F and G functions have already been dis- 


cussed, 


Now that we have constructively arrived at a general expression for sur- 
faces that have a prescribed boundary vector and a prescribed boundary slope 
vector, it might be interesting to apply a proof to a conjectured higher order 
surface equation in which not only boundaries, boundary slopes, but also 

. boundary second derivatives are vector quantities under control. 


We postulate, therefore, that by analogy the surface equation is 


uw =- (-2 Fou Fyu Gu Gu Hou Hu] 


0 1 


1 


It represents a surface patch whose vectors of coordinates, slope, and curva- 
ture as well, are everywhere arbitrary along its boundaries, The first column 
and first row of the 7 x 7 matrix represent these boundary conditions; the re- 
mainder of the matrix is redundant, since the quantities this partition contains 


must all come from the column and row by differentiation. 


We can test this equation by seeing whether it contains a boundary curve, 
To this end, set u = 0, so that we check whether it contains the boundary (0w). 
We obtain, invoking the stipulations on the F, G, and H functions, 
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Ke 1 


In the boundary matrix we have omitted irrelevant terms, because of the 
zero's in the pre-multiplying vector. We obtain, by performing the multipli- 


cation, 


= Ow, which is the hoped-for identity. 


We can next try to see whether the equation also conforms to the boundary 
second derivative conditions. It will be more convenient in what follows to 


introduce some abbreviated notation, 
Set 
[tu] ~ [-1 Fou Fu Gu Gu Hou Hl 


and a similar expression for [tw | é 
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Set the 7 x 7 boundary condition matrix equal to [3] . 
With these abbreviations, the surface equation is 
ow e- [a] fl foo] 7. 
We differentiate with respect to u: 
we (fed f+ Es} fe)) fo] 7 
and again; 
mae (fem) [5] +2 [ey] fa] + [8] fall) 7 
We wish to investigate the right hand side of this equation for u = 0, thatis, 
for UW = Owia’ The blending function vector and its derivatives become 


mq] = [0 0 0 001 g . 
f' 1 0 0 0 
[to] = [-1 10000 0 


As for the first and second partial derivatives of the [a] matrix, all 
elements of [B | and [B | vanish except for those in the top row. 
u uu 


Then 


N 
— 

we 

So 
poner) 
—s 


B,| = [ the null vector] 
[ro] [s., -[ 0 ORs Or. ~00 -01 ~-00 -01 


The sum of these vectors is evidently 
Finally , 


= OW a , a8 expected, 
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We have shown that the extended surface equation satisfies the second 
derivative boundary conditions. In a similar way it can be shown to satisfy the 
first derivative boundary conditions, but this was skipped in favor of the proof 
for the higher derivative, since the procedure exhibits a few interesting points. 


By analogy we could construct matrix products to represent surfaces 
which satisfy even higher derivative conditions across boundaries. 


3.5 BOUNDARY CURVES 


It is often convenient to use particular boundary curve functions defined 
by the curve end—points and end-point tangent vectors. We can use the blending 
functions themselves to define such curves. For example, the u0 boundary 


curve can be described by the equation 


ud = [Fu Fu Gou cu] 00 


where the column vector contains the end-point information. We observe that 

the row vector becomes (a 0 0 o| ‘when u = 0; it becomes [o 1 0 o| F 
whenu=1. Again, if we take derivatives of this row vector with respect to u 

we obtain [F)'u Fu Gu G,'u] and this becomes [0 0 1 0] foru=0, 
and it becomes [0 0 0 1] foru=1, 


With this behavior of the row vector, it is easy to see that the equation 


does indeed represent a curve satisfying the end-point conditions. 


The matrix form of the surface equation has been shown to be 


iw fa Fu Fyu Gou Gu] 0 uo oul ud ul, |[-1 
Ow 00 01 00 01 Fw 

w w 0 
Ow 00 01 00 OL G.w 

u uw uw. 0 
lw 10 11 10 11 G.w 
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Now when, in computing uw, we perform the matrix multiplication from 


the left. we have for the second column 


[-1 Fou Pyu Gu Gul Tuo 
“00 
10 
| 00 
| u 
; 10% 
> yet 
= -u0 - [F u Flu Gu G ul "90 | = ul - ud = 0, 
Q 1 0 li | 
10 ! 
00 
u 
10 
u 


Now if similarly ul, HO... and ul, are functions of the same kind, their 
corresponding column products vanish just as in the case of uO. Accordingly, 


the resulting product of the three matrices has the form 


uw ks 0 0 Of “gh ap, 
J 


ner wi 


) 
Ce 


where P is the product of the row vector and the first column of the matrix, or 


, 2 s 3 cl r 
P= Ee Fou Pu Gou Gu; lw | 


0 


a [F u cr Gye Gul Ow 
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If again the elements of the column vector, Ow, lw, ow and Iw, are 
likewise functions described as outlined, we can write typically 


Ow = [00 01 00 01.) Fiw 
w Ww 


Ow = 00 01 00 01 Fow 
lw 10 11 10 11 Fw 
Ow 00 01 00 01 G,w 
u u u uw uw 0 
lw 10 11 10 ll Giw 
u u u uw uw 1 


When we substitute this result for the column vector, we obtain the surface 


equation 
uw = [ Fu Fu Gou Gu] Fow 
Fw 
Gow 
G.w 


This is a particularly convenient form for computation. The 4x 4 
matrix contains nothing but information about the corner coordinates, corner 
slopes, and corner twists; all entries are constants, and the partitions of the 
matrix systematically group these quantities. The leading row vector and 
the trailing column vector are transposes of one another, (but with different 


arguments, of course.) 


We shall refer to the 4 x 4 matrix as the "boundary condition" matrix, 
and shall assign to it the symbol B, so that the matric equation for the surface 
could be written 
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uw = [ Fou Fu Gou Gu] B Fow 


1 


Gow 
G.w 
It must be remembered that each of the entries in B is a three-vector, 


whose components are x, y, and z coordinates and slopes and twists. This 


means that B is really a tensor. 
3.6 BLENDING FUNCTIONS 


We can relate the blending function vector to a so-called basis vector 


in the following way. Let [u, u ] be a vector whose elements are a set 


243 "4 
of linearly independent functions of the variable u. Then we can postulate 


the existence of a matrix M such that 


[ Fou Fu Gou Gu] = fu, Uy Us U4] M. 


To evaluate the M matrix, we substitute u = 0, u = 1 on both left and right 
hand sides of the equation. Then we take derivatives of both sides, and 


again substitute u = 0 andu = 1. There results 


F)0 F,0 Gy? G,0 = o 0, 04 0, 
Fol Fil Gl Gl 1, 1, 1, 1, - 
t t t t Ul + ' ' 
F? Fe Ge G0 0 0 0 o 
i 1 1 1 t £ ' ' 
Fol Fil Gol Gl 1, 1, 1, 7 
The matrix on the left is the identity matrix 1000 by virtue 
0 1 0 0 
001 0 
0001 
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of the stipulations on the blending functions and their derivatives. From this, 
we conclude that 0 00 0 


-1 


and we need only find the inverse of this matrix (if possible) to obtain M. 


(in the matrix, the notation 1, means d(u,) typically. ) 


du u=1 


In the next section we shall for the first time be specific about the 
Uy Ug Uy ], but it is interesting and important to realize 
that so far in the discussion nothing has been said to diminish the generality 
of the mathematical structure. It is hoped that the reader will not lose 
sight of the fact that the surface equations in their several forms can be 


implemented in many ways. We propose to develop one such implementation ~ 


basis vector [ a 


in detail, but it is only one of many. 
3.7 CUBIC BASIS VECTOR 
Let the basis vector be 
2 3 4 


fu, 4 u ep Se ae at), 


The vector on the right contains four specifically chosen linearly independent 
functions of u, the powers of u, and when multiplied by a coefficient vector 


yields cubic polynomials: 


tue wu lj = Au? + Bu’ + Cu 4D. 


59am Pp 
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By the reasoning of the last section, we have for this vector basis, 


ot 


KH » & © 
a — 


0 
1 
0 
2 


ao oS = © 


from which we can obtain the desired inverse 


22 1 21 
M= {-3 3 -2 =-1 
0 0 1 0 
1 0 0 90 


Now we can write 


3.2 
[Fu Fyu Gou Gu] =fueu ou ij M. 


We shall abbreviate the notation for the basis vector in what follows. We 
shall write 


tue wou ips U 


and tw? w wlj.= W. 
The matrice surface equation 


uw = [ Fou Fu Gou Gu) B Fow 


0 
Fw 
Gow 
G,w 
now becomes, simply and compactly, 
uw = UMB mt wi, (Superscript t means transpose. ) 


If U and W are cubic basis vectors, then the surface patch is the so-called 
bi-cubic surface. Such surfaces are very easy to compute, particularly 
since the basis vector is so easy to evaluate. In passing it is important to 
remark that the above compact surface equation is not limited to cubics; 
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U and W are not restricted to cubic basis vectors, and M is simply the matrix 
that generates the appropriate set of blending functions. Among other pos- 
sibilities, U and W might be higher order polynomial basis vectors; or they 
might be any set of linearly independent functions. Provided the associated 


ut matrix has an inverse, these basis vectors are acceptable. 


We can write, for w held fixed, an expression for a u ~ varying curve 
on the surface: 
t .t 
uw = U(MBM W) = UA 


where A is a column vector of constant coefficients. We can write a simi~ 
lar expression for u held fixed and w varying. The matrix product mpm 

is the same in either case. This suggests that for any surface patch this 
product should be evaluated first; thereafter, we can either obtain u-varying, 
w-constant curves or w-varying, u-constant curves in an obvious way. 


We shall investigate another basis vector that is composed of another 
set of linearly independent functions (not powers of u) in a later article. 


3.8 DIFFERENCE EQUATIONS 


If the basis vectors are polynomial bases, we can invoke the techniques 
of finite differences to calculate points on the surface patch. 


Consider the matrix 


Le= 1000 
110 0 
1141 =(0 
1111 
Then a a 
L b = a+b 
e atbe+c 
d atb+c+d 


If a, b, c, d are respectively third, second, first and zero-order differences 


of the cubic n°, then the column matrix on the right of the equation represents 


os 
3) 
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the corresponding differences for the cubic (n + »’. The differences for 


(n+ Ky? are given by 


1 
c e where is means k successive 
b 


Q 


Q 


i 
_d 
multiplications by the L matrix. 


When n = 0, we can easily find that 


for cubics. 


Using this, 


cs [ 0001 ] i ia 
=6: 2 


In this expression, the vector [0 0 0 1] serves to select the bottom element 


of the resulting column vector after k multiplications by L. 


By extension, we can write the more complete statement 


3 k 


(ee bape few aye 6 0 0 o| 
= ee o | 
de tcic oe i 
0 0 0 i | 


We shall call the [0 0 0 1} vector 1; and we shall call the 4 x 4 matrix N, 


so that 


2 
ek? ea = ra 
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Now let the usual parametric variable u be represented by u = k5, where 
k = 0,1, 2... .. and where 8 is an increment size. Then 


tue wou i] = rk? i k 1] §3 
§2 


Call this last square matrix 4. Then 
Us irk NA. 


This expression states that we may step along the u parameter, in 3 incre- 
ments, by successive multiplications by the L matrix, and thus evaluate 
the U vector at these steps. In order for u to go from 0 to 1, k must start 
at 0 and go to —: since k 5 = 1, 


We can also write, for the W parameter 


We ik wa 


The surface equation in difference form can now be written out in 
full: 


uw = inna MB mt A nt i re 


Call the partial product NA M B Mm’ At nt = Soo # square 4x 4 matrix. 


Then L S50 = S10 @ new square matrix, 


and Us 


00 = 85° after j multiplications. 


We remember that for any column of Soo the multiplication by L is a process 


of cumulative addition, as shown by 


& 
>) 


al 
o 
tt 
fo 
+ 
o 


io%) 
~ 
+ 
log 
+ 
Q 
+ 
a. 
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We can write in general that 


t 
L L = 
Sy Bed, kel 


where the new square matrix is obtained from the old by cumulative addition 

of column elements, followed by cumulative addition of row elements. These 
operations are furthermore commutative, which means that we obtain the same 
result if we first add row elements and then afterward add column elements: 


We have, finally, that atu= 0, w=0 
the surface equation is 


00 = 189 1 * and in general 
uw = 18), 1° where u = jd, w=kd, and S,. has been formed 


from S 00 by j column additions and k row additions. 


This obviously furnishes an extremely simple way to generate discrete points 
on a surface patch. The pre-multiplier 1 has the effect of selecting the bottom 
row of Si and similarly the post-multiplier i has the effect of selecting the 


last column of Se The bottom right hand corner element of Six is the value 


of the coordinate for a point on the surface, atu=j5, w = k5. 

Consider 1 s; 0° This represents the row vector obtained after j cumula- 
tive addition operations have been performed on the columns of Soo: The 
right hand element of this vector is the value of the coordinate atu = j9, 
w=0. We can hold u fixed and step out successive values of the coordinate 
for w - varying, simply by cumulative addition on this row vector alone. In 


this case, the resulting right hand element is the marching coordinate value. 


An analogous remark can be made for the product Sox 1 af This is a 
column vector, and successive cumulative additions of its elements marches 


out values of the surface coordinate for w = k® fixed, and u varying. 


Although the arithmetic of the foregoing difference method is very 
attractive, it possesses certain drawbacks that must be made explicit. The 
coordinate values are precise if and only if no truncation error whatever is 


allowed in the arithmetic. Error is cumulative, and the least departure from 
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a precise number at the start will rapidly propagate. In the Amatrix there 
are numbers of 8 magnitude; for 100 calculated points along a curve, this 
calls for 6 decimal digits, all exact, or 18 binary bits. For display purposes, 
100 points along a curve are adequate, but for engineering purposes this 
information is too sparse. To calculate intermediate points, the 4 matrices 
need to be changed and the S matrix recalculated. Furthermore, it is clear 
that we soon reach an upper limit on the number of available bits in the 
computer word, becasue of the rapid growth of 83 


We can of course use the difference technique to calculate coefficients 
for a cubic and then calculate points using them. For consider 


uw = UMBM- at x? I, 


In this, the partial product M B ma at nt ur represents a square 


matrix whose last column consists of the coefficients for the cubic at fixed 
w = k8, with u varying: If S, is this square matrix, 


w= US 175 U A 
B 
Cc 
D 


The cubic can be generated by digital methods for any u, or it could be 
generated quite easily by analog differential analyses hardware. The inte- 
grators of the analog differential analyses are loaded from the values of the 
column coefficient vector, and it is then unclamped and allowed to generate 
the curve. Meanwhile the digital machine can perform a cumulative addition 
on the rows of §. to yield Ss, ‘o The last column of this new matrix then 
contains the new coefficients for the curve of w = (k +1)5, ready to be 
loaded into the integrators. 


Recent developments in hardware, particularly a hybrid digital-analog 
multiplying device, may make it possible to generate surfaces for display 
directly from the indicial forms of the equation, or from the matrix form 

uw=UMBM wY 
without the necessity for recourse to this last difference technique. It is 
however interesting to know that the method exists. , 
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SECTION IV 
HY PERSURFACES — HIGHER DIMENSIONS 


We can readily extend the surface equation to describe hyper-surfaces 
immersed in hyper-space. For this purpose we shall introduce a slight variant 
on our notation. We shall write typically 


u 
i 


Fu. i=Oorl, 


This will be a standard replacement for the blending function notation. The 


stipulations on the Fi are as before, so that ifu =a, a =0or1, we can write 


a 
i 


it 


0 when a #i 


a 
i 


1 when a =i. 


For slope continuity across boundaries, typically 
t 
- = 0 where this symbol means the first derivative of the biending 

function with the argument = a. 


For higher order continuity across boundaries, the additional stipulations 
on the blending fucntions are the same as for ordinary surfaces with two degrees 


of freedom and have already been discussed. 


The general surface equation for hyper-space is, in indicial form, 


“(uvw...) =k. TR: : 
+ (ivk ie 
+ (ijw . iy : 


- (Nel) (ijk. . eae a 
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In this equation, N is the number of independent parameters in (uv w...); it is 
the number of degrees of freedom of a point on the hyper-surface. The indices 


i, j, k etc. can take on only the values 0 or 1. 


Let us proceed to prove that this surface contains a boundary, say for 


example the boundary (u00.. .). We hope that the following equality holds: 


. ; 00 
quod...) SARS ink sae 8 
. u 0 
+(GU0k.. Vike 
uO 
# (kJ Oee tars .e 
Gj a 
+ 
A rn ud 0 
CASA MI A SE Sg ie Paes 


The last term in this expression is non-vanishing if and only if all indices other 


than i are zero, i.e., j= 0, k =0, etc. We can accordingly rewrite this term 


as 
(N-1) G@00...),. 


Next consider the second term on the right: 


7 u 0 
qgGok.. D5 ko? 
It is non-vanishing if and only if k = 0, etc. 


We can accordingly rewrite it as 


(i00...),. 


A Similar consideration applies to 


ud 


et which also becomes 


(ijO...) 


GOO; e045: 
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There are evidently N-1 such terms, all identical, and they are removed by the 
last term. All that is left is the term 


but since j = 0 and k = 0, we finally have (u0 0...) onthe right. This establishes 
the identify, and the surface equation has thus been shown to contain this boundary 
curve. It is trivial to show that the surface contains all boundaries, and is defined 
by them. 


We can also show that the hyper-surface contains boundary surfaces 
of lower order. We shall content ourselves with the case for N = 3, and show 
that it contains surfaces for N = 2 which are identical with our ordinary surfaces. 


We have 


v 


(uv w) =WiW yy 


+(vkyiy 


+(Chiwy 


-2d WEY 


Set v=0. Then substituting, and retaining only non-vanishing terms, 


(which means that } must be replaced by 0 whenever it occurs, and j = : = 1) 
(uOw)=(u0)” 


+(i0k) Sy 
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or 


(ow =(U0k) ” 


* (Ow); 


3 uw 
- (i 0k) ik: 
This is the two-degree-of-freedom surface 


Ww u uw 
= (uk + (iw -(i 
(uw) (uk) k (iw) i (ik) | k 
We shall next consider the slope vector of such a hyper-surface. 


partial derivatives with respect to one of the variables, say u, and get 


ahs vw 
CURING ace Di = COURSE eas ae 


cet at ul w 
ee ae ea 
se pace. ul v 
(ijw Maa 
ulvw 
~ (N- ijk...), Sw 
(N-1) (ijk ae 
In this, set u = 0, 
vw 
foe Gale's =(Ojk... Wh ceca 
vw...) = (05 ee 
O'w 
*(ivk.e.w. 
ey vik 
O'v 


PL Wes fe Gaus kot 


-(N-l) (ijk...). 


We take 
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or 


(Ovw.. dy 7@ik. Z Het: Sato 


since all other terms vanish by virtue of : =0. 


This result is analogous to the one obtained for boundaries of ordinary 
surfaces; it says that the slope anywhere on a boundary is a function only of the 
slopes at the "ends" of the boundary, and are otherwise independent of the 
boundary shapes, Slope continuity across boundaries is a consequence. 


The hyper-surface equation just developed is defined by ordinary curves, 
or single-degree-of freedom boundaries; we can also write a hyper-surface 
equation for N degrees of freedom, defined by boundaries with N-1 degrees of 
freedom. We shall exhibit the result for N = 3: 


(u v w) =(ivw)/ 
*wiw, 
+uvky) 


ato]. 


uw 
-@vkyiy 


The proof that this space contains, for example, the boundary subspace 


(0 v w) follows the preceding proofs in principle and will not be carried out. 
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SECTION V 
SURFACE NORMAL VECTORS 


5.1 GENERAL SLOPE CONTINUITY CRITERIA 


The surface normal vector furnishes a convenient mechanism for the in- 
vestigation of general criteria for continuity of surface slope across boundaries 
between surface patches. I will be seen that the continuity conditions already 
established are much stronger than are necessary, but that they are expedient. 


Put \ 


ue [x94] 


w= [ay xy 
for the tangent vectors of a surface patch at some point. Let us assume that 
another adjacent surface has a common boundary curve along u = constant, w 


varying, so that Wis common to both patches. Let the parameter for this second 
patch be v, and for its tangent vector put 


Ve lie yy z | . 


COMMON 
BOUNDARY 
CURVE 
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The two surfaces will be continuous in slope across the boundary at the point in 


question in case the three vectors U, W, and V are coplanar there. 
The surface normal vector for one surface is 
N=UxW. 
If V is perpendicular to N, 
V-N=0, 


But then in this case, the three vectors U, W and V are coplanar since they are 


all perpendicular to N. 
In detail, this gives 


N=UxW= [5,4, J, ] 


a vector whose components are the familiar Jacobians, and 


V°N = vue 
: [x, Vy z, an = x, Vy 2, a 
a x, ye zy 
a wy Vw tw 


(This is the so-called "scalar triple product" of the vectors.) Thus the vanishing 
of the determinant of the matrix of the three tangent vectors is the general condi- 
tion for slope continuity between two patches, at any point on their common 
boundary. 

This equation also shows that we may have slope continuity of surfaces 
even though the curvilinear coordinates of the two surfaces are not slope - 


continuous across the boundary. 


If the tangent vectors U and V are equal everywhere along the boundary 
curve, the determinant is sure to vanish; similarly if the tangent vectors U and 
V are scalar multiples of one another, even when the scalar multiplier is a 


variable quantity. 
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GENERAL CONSTRUCTION - TANGENT SURFACES 


Suppose that a surface A already exists, defined by the parametric vector 


equation 


A= [x 6) vib 9) 26 @)-] 
Let a be a curve on the surface; it is always possible to write the vector equation 


for a in either of two forms: 


ap )= (x) vib) xo) 


or 
a @)= [x 6) a) 2 8)] 


Suppose we wish to attach a surface B to surface A, in such a way as to make 
curve a common to both surfaces, and suppose furthermore that we wish to 


maintain slope continuity across this mutual boundary. 


We shall consider curve a to be the boundary (0w) of the B surface. We 
are at liberty to design, arbitrarily, a projection of the other three boundaries, 
(uO), (ul), and (lw), Say for example that we design these curves in the xy pro- 
jection. Then the curves represent the x and y components of their coordinate 


vectors. 
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We are now ready to obtain the missing z component of the tangent vector 
across (Ow). We first compute the surface normal to A along curve a. For this 


purpose we can use any one of the expressions 
PE), x (h OQ 
or Oa), x O)g 


or ab), x (P8)g 


Each expression yields a surface normal N; the three results are identical, 
We can evaluate this surface normal vector at any point on (Ow) since we have a 


correspondence between w and the variables ¢ and @ . 
We also have the equation 


(Ow), ni = 0. This is the familiar condition for surface slope 


continuity. 
Let N= [a b e] after evaluation at w. 
Then the equation becomes 
[xiow), y(ow), 2(0w), ] [a 
b =0 
c 
We already have the x and y components of this equation, and can solve for the 
z component: 


z(0w), = itachi a otal 


“Cc 
This z component has a magnitude that ensures that the complete vector 
(Ow), is coplanar with surface A at w. Hence (Ow), is the desired tangent vector 


of surface B across (0w). 
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We next find the z tangent vector components 2(00) and z(01), from (Ow) 
and use them in the equation 


2(0w) + = 2(00) Fow + 2(01) F,w. 


Here the 1 superscript indicates that this is an intermediate result; it is the 
intrinsic boundary tangent vector for the F-type surface, and does not yet match 
the 2(0W) vector function obtained from the A surface. 


Accordingly, we must add to the F-type surface a G~type tangent vector 
correction surface, so as to make the combination have the desired slope along 


(Ow). 


This G-type correction surface is, as we have already shown, 


c c c 
{uw)” = [1G,u Gu] 0 uo, uly 1 
ow’  =00° -01° G.w 
u uw uw 0 
lwo =10" “1° | jaw! . 
uw uw 1 


The superscript c indicates that this is a correction surface. 


Slope correction is necessary only along the boundary (Ow); we can enter 


the value for owe in the matrix, but the other entries must be looked at in detail. 


We have, for the slope correction across (Ow), 


ow’ = ow - ow! , 
u u u 
These latter two quantities have already been found for the z component, and so 


c, 
a is known, 


Consider lwe . This is at a free boundary, (lw), remote from (Ow), and 
we can set it equal to 0. Then 10° and 11 are both zero also. On the other 
hand wo" and uly are connected to (Ow) at (00) and (01) and so we must specify 
them in such a way as to satisfy the conditions at these points. Elsewhere, they 


too are arbitrary. 
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We write the G function expression: 


(uw) = [1 Gyu Gu] 0 uo, ul, 1 
Ow -00 ~01 G,w 

u uw uw 0 
0 0 0 Giw 


(we omit the c super-script temporarily. ) 


Performing the first multiplication, we have 
1 


ees [ow, Gu | ca “0 or | La 01, SoU] G.w 


0 
G w 


Consider the element ud - 00 woot in the row vector. Since ud is 
arbitrary, it would be convenient to choose it so as to make the entire vector 
element vanish. We therefore write 


u0Q = 00 Gu, 
w uw 


0 


Then 00. 0 as it should, and 


’ 
u, 


since ud =00 G 
uw 0 


wu 


00 = 00 _ as it should. 
wu uw 


Similarly, we may set 


ul = 01. Gu. 
Ww uw 0 


The result of these choices of wd, and us is to reduce the G equation to 


uw’ = Ow’ Gu. 
u 0 


This represents the correction surface z component which must be added to the 
z component of the F surface in order to obtain slope continuity along (Ow) between 


the given A surface and the designed B surface. 
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5.2 ADJACENT-PATCH SLOPE CONTINUITY 


It is sometimes desirable to define the boundary curves for two adjacent 
patches so that at the junction between the curves the tangent vectors have the 
same direction but are of different magnitudes. This is particularly useful when 
the boundary curves are parametric cubics, because then the magnitudes of the 
tangent-vectors at the end points control the behavior of the curve segment. 


ul 


Ow 


do ud 
00 


As a specific case, consider the boundary (Ow) common to two patches; let 
the tangent vectors at (00) and (01) for the first patch be 004 and 01 p? and let 
the tangent vectors for the next patch be 00. and 01, 


If the tangent vectors have the same direction, they are scalar multiples of one 


another, 
or 00. =m Ss 
and 01, =n Oly ‘ 


Suppose that the patch @ w) already exists. We need to obtain an appropriate 

expression for (uw) so as to match surface slopes across (Ow). By the results 
of the preceeding article, we can accomplish this in very general ways, but in 
our present case let us make a special requirement on the tangent vectors: let 
us assume that everywhere across (Ow) the tangent vectors have the same 


direction, and differ only in magnitude. Then for any w, 


ow, = r Ow 
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where \ is ascalar. We know that \ takes the value m at w=0, and takes the 
value n atw=1, and we conclude that )\ therefore must be a scalar function 


of w. 


With this relationship between (ow ) and (ow) the vector cross product is 
always the null vector: 


ow), x Ow = [0 0 0] 
Hence the scalar triple product of (Ow), (Ow) f° and (Ow) vanishes for any 
(Ow). This ensures that the two surfaces will be continuous in slope across 


(Ow) for any shape of (Ow) and for any » = > (w) that has the proper behavior 


at w=0 and w=1, 
We could for example take 
\(w) = m(l-w) + nw. 


This is a linear variation of \ with respect to w. It has one disadvantage, 
however, in that it introduces un-wanted cross derivatives or twists at (00) 


and (01). In order to avoid this, we might use 


d»\(w) = m Fow + n Fw. 


Then the required slope function across (Ow) for the (uw) patch is 


ow = fm Fw + n Fw) Or: 


We can check to find the cross derivatives introduced by this relationship. 
The cross-derivative is obtained by differentiating with respect to w, and 
yields 


u 


t ' 
Ow (m Fyw + a Fiw) ows 


uw 0 


+ mM Fow + n Fw) Owed 


At (00), 00 =m 00 
uw 


pw 


and at (01), 01 = n Ow 
uw 


ow 


54 SECTION V 


This shows that the \ function does not introduce additional twists at the 
corners of the patch, beyond, of course, those already inherent in the (@w) 


surface. 


If the (uw) surface already exists, defined perhaps by an F-type equation, 
1 
its intrinsic tangent vector across (Ow) is a known function of (Ow), say ow 
where the 1 superscript indicates that it is an intermediate result. Then, 


as before, the correction of slope is 
Ow” = ow, - Owl. The correction surface 
is (uw) = Ow Gu. 


When this correction surface is added to the original surface, the combination 
will be continuous in slope with the (#w) surface across (Ow). The u and 
curvilinear coordinates of the two patches will be continuous in slope across 


(Ow), but their tangent vectors will be different in magnitude. 
5.3 APPLICATIONS 
Let dU represent a differential vector, so that 


au = [ dx dy dz}, in which 


_ 2x ox 
dx = 77 du + oy oY 
ay = SX au + 3X aw, 


2 22 ay 4 OF 
dz dy et 5 Ww. 


If dU is tangent to a u-varying, w-fixed curve, these become 


dx = Tq =x du 
du u 
dy = a du. 
dz = zo du, since dw = 0. 


Hence dU = [x y 2 ] du. 
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Similarly, if dW is a differential vector tangent to a w-varying, u-fixed 


curve, 


dw = [x % 2y] dw. 


The normal differential vector at a point of the surface will be given by 
the vector cross product: 


dn = dU x dw 
Yuu 24 u *%u 
= du dw. 
Vw aw aww Ww 


The determinants that comprise:the elements of the vector are the 


jabobians Je de and J, so that we may write 
dn = [3.3.3] du aw. 
x y Z 


The magnitude of dn is equal to the differential area of the elemental parallelo- 
gram described by dU and dW. This magnitude is 


| an | =afdn:dn r./ [ac] [an] T) 
du dw VIi2 + Jy2 + v2 . 


From this, we can construct an algorithm for finding surface areas of patches; 


we simply perform numerical integration of the expression 
‘1 1 
ieee ging ~/J2+359 +d dudw. 
0 0 ‘ y 2 


Again, if N is the unit normal vector to the surface at a point, then 


dn 


N | dn | , from which 


where 8 


u 


J + : 
J2 Jy2 rag 
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The quantities 5, Jy, Jn are the coordinate components a, b and c 
s S s 


of the surface normal. 


If the surface is to be manufactured by milling with a ball-end cutter of 
radius R, then the cutter-center vector [x, v. zy ] is related to the surface 
vector [x y z| by the simple expression 


[x, oe z,] = [x y z| +R [ab 7 
This describes a "parallel" surface spaced a distance R away from the 


designed surface. 


The normal vector can also be used to calculate volumes enclosed by 


surface patches and planes, as follows. 


Suppose we wish to calculate the volume contained between a surface 


patch and the xz plane. 


XZ plane 


We can imagine the volume broken up into a number of slender prisms 
whose axes are all parallel to the y axis (and perpendicular to the xz plane.) 
The area of the base of one of these prisms is the projection of the small 


element of surface area, or 


dA = J_ du dw. 
y y 
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The volume of this prism is 


qdVo= ydA = aoe du dw 


A f 
vVoa] ! yd. du dw. 
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SECTION VI 
CORNER TWIST VECTORS 


6.1 THE QUASI-SPHERE 


It is possible to choose a parametric cubic that very nearly approximates 
a circle for one quadrant. We shall go into detail about this shortly; intuition 
suggests that similarly we ought to be able to construct an approximation of an 
octant of a sphere by means of a bicubic surface, bounded by these approxi- 


mations to circles. 


For the circle approximation, let us assume that we will be content to 
make the quasi-circle pass through a point on the true circle atu = 1/2. (This 
is not the best possible approximation, but it yields quite good results and the 


arithmetic is simple.) 


We shall assume a circle of unit 
radius, centered at the origin, with 
end-point values of the parameter u 
as shown. The tangent vectors are 
symmetric, but have yet undefined 


magnitudes. 


We have 


x= [u uiou i| M {| 0 Where the column vector on the 


right represents the end conditions 


- 


for the curve. 
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When u = 1/2, this becomes 


noe [1 2 a8] fo  -2 l 1 0 
-3 a l 
0 0 l 0 a 
| 0 0 0 0 
> [ ] | 
ee (li Bae a -Q+a 
8 
3 - 2a 
a 
L o 


i 
| 


= (-2+a+6-4a+4a) = > +a) 


a= 8x-4. Butatu=1/2, x= 2. since it is a point on the circle 


(by symmetry, at Zi 


Hence a = 4 (2 - 1). This is the required magnitude of the two tangent 
vectors at u=0 and u = 1. (Calculation reveals that the quasi-circle has a 


radius of about 1.00016 at u = 1/3 (at? or 30°) so it is a good approximation.) 


We now establish a coordinate system for the sphere, and show its 


boundary curves. 


10. 1 


ud = degenerate curve 
boundary 
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The boundary curves Ow, lw, and ul are all unit circles; the boundary u0 


however is a degenerate circle, and appears as a point. 


We shall first investigate the z component of the uw surface vector 


- er 
= : 1 1 
z(uw) [Fou Fu Gu Gul 00 0 00, 01, F yw 
10 11 10 ll Flw 
Ww Ww 1 
00 01 00 01 Gwi] 
u u uw uw 0 | 
10 ll 10 li Gw! 
L u u uw uw] 1 3 
= 4 -a | ; 
[Fu Fu Gou Gu] 0 | o a Fowl. 
1 0 0 -a F_w 
0 0 0) 0} Gow 
0 0 0 0 G.w 


When we perform the first multiplication, we obtain 


z(uw) = [Foe ~ Fou | 0 | 0 | -a (Fu + Fw] Fw 


But Fuu + Fyu = 1, by virtue of the definition of the F functions. 


Hence 


ats cae 
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The curvilinear coordinates for w constant thus yield constant z; this implies 
that z is independent of the other variable u. and the w curves are plane curves. 


They must of course be quasi-circles. 


We have obtained the value of the number a by investigating a unit circle. 
For a circle of radius r, the tangent vector magnitudes must be proportional to 
r, or equal to ra. We can find these radii for various values of w from either 


the quasi-circle y(Ow) or x(lw). 


We have 


x(1lw) = [10 11 (10 PT ] Fow 
W wi 
Fow 
1 
Giw 
0 
Giw 
and 
2 re dé 
vow) [00 01 00 ot. | Fw 
Ww W. 0 
Flw 
1 V 
Gow 
0 
G 
es au 
in either case, 
r(Ow) = fo Loa 0 | [i w| 
‘ J 0 
Fow 
] 
Gow 
G 
i“ 
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For w fixed, the x and y coordinates of a quasi-circle are given by 


x = [Fu Fu Gyu Gul 0 


r 


and 


q EF 
[Fou 1° ihe Gu] Yr 


ke 
tit 


where r is a function of w, shown above. 


But 
cr - r 
0 0 0 0 0] row 
} iT 
r 0 1 a 0 Ps i 
= : 
ar 0 a a 0 j Gy | 
| 0 0 0 0 0 ee 
and 


[ Yr [9 1 a 0 o 
| 0 | 0 Q 0 0 FW 
0 0 0 ) Gow 
| | 0 0 | 6 
2 | 
| -ar| 0 -a -a 0 | |G w 


We obtain these last results by simply writing rows in the 4x 4 matrix that 
correspond to elements in the vector of the left side. When we combine 


results, we have 
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x(uw) = [Foe Fu Guu Gu] 


yuw) = [Fou Fu Gu Gu] 


The equation for z(uw) has already been shown. 


The striking thing about the B matrix as it appears in these equations is 
that is has non-zero entries in the bottom right partition, 


00 ol 
uw uw 
10 11 By comparison, we see 
uw uw 
that x(00 = a 
uw) 
and (10 = <a 
yO.) = -a. 


These are the cross-derivatives at the corners 00 and 10. All other cross 
derivatives vanish. We shall refer to these cross-derivatives as "twists" of 
the surface; UW aw is the twist vector at a generalized point on the surface. 


6.2 THE EXACT SPHERE 


A bi-cubic surface cannot fit a sphere exactly, and it would be interesting 
to see whether by an appropriate choice of F and G functions other than cubics, 


an exact equation can be constructed. 
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The functions 


and 


F.0 = Fy = 0 
' ' 
Foo = Fol = 0 
as well as 
FO = 1 Fil = 0 
F,0 0 Fil ane 


We shall choose these functions, and determine appropriate G functions 


so that the equation 


(uy) = [Foe Pou cot Gu] 


represents an exact circle, and not the approximation of the last section, 


We already have the well-known parametric equations for a circle; 
x ~ sin au 


y = cos au where a 


mf 
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X 
We can compare the x equation with 
2 2 
x(uy - (cos au sin au Gu Gul 0] 
1 
a 
0 
2 
= sin au - aGu 
0 
2 
Then sinau = sin au + GG us whence 
C 1 2 
7,U +» ~ (sinau - sin au). 
49 aie ) 
Again, 
2 2 | 7 
y = cosau [ cos au sin au | Guu Guy fe 
, 0 
Lo 
| 
+a 
2 
cos au = cos au - aGu 
whence 


iL 2 
Glu = a (cos au - cosau). 
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We can easily verify that these G functions satistv the same stipulations 
as the cubic G functions: 


3,0 = Gd = ( 
Gy 0 


= = = 0 
GO Gl 0 


' 
30-1 = 
G9 Gale 


a 1 


G0 = 0 Gl = <1" 


If we now use the same boundary value matrices as were used in the 


: . , : Tw ae 
previous case of the quasi-sphere, but with a = ae throughout, we obtain for 


the z component of the surface vector 


z(uw) Eat -2 Gyw ; 


This is. with the new F and G functions 


2 2 
z(uw) ~ cos” aw - (cos” aw - cos aw) 


cos aw. 


As before. this shows that the z coordinate of the surface is independent 
of u: the w curves are plane curves. and they are indeed circles. Their radii 
are given by 


x(lw) lio Tl 610 1] Fy. 


oo 
= } 


= Fw - ae 


2 2 
sin aw - (sinaw - sin aw) 


or r = sinaw., 
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Then 
x(uw) = [Fu Fu Gou Gu] 0 


ar 


= r(Fju + aGou) = sin aw sin au, 


and 


y(uw) [rye Fou Gyu Gul r 


-ar 


= r(Fou - aG,u) 


2 2 
= sin aw (cos au - cos au ~- cos au) 


— sin aw cos au. 


The resulting parametric equations, when collected, are 


x = Sinaw sin au 
y = sin aw cos au 
Z = cos aw 


and these are well-known. 


This demonstrates that the sphere is a special case of the general surface 


equation, provided the blending functions are suitably chosen. 


The F and G functions are by inspection, seen to be linear combinations of 


the linearly independent functions of u, 


2 an F 
[cos au sin au cos au sin au] a 


and this may be taken as an appropriate basis vector. Then 


2 2 
[ Fou Fyu Sy Gu = [cos au sin au cosau sinau | [a] 
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where the M matrix is, in this case, 


1 0 0 


[ 1 
a 
0 i ose See 
a 
M = 1 
0 0 o - 
a 
l 
0 0 = 0 


Incidentally, its inverse is 


1 0 r 0 
eis 0 1 0 1 
: aan ne) 0 0 a 

0 0 -a * 


This empty page was substituted for a 
blank page in the original document. 


71 


SECTION VI 
RULED SURFACES 


The locus of straight lines connecting corresponding points on two curves 


(uO) and (ul) is a surface. The lines are called "rulings" of the surface. 


ul 


ud 


The equation for a generalized line of the surface is also the equation for 
the surface: 
(uw) = (uljw - (u0)w + (ud). 
This is equivalent to 
(uw) = (uljw + (ud) (1-w). 
From this equation we obtain the derivatives 


(uw), = (ul) w + (uO) -w) 
(uy = (wl) - (ud) 


(uw), = (ul), - (ud), . 


7.1 DEVELOPABLE SURFACES 


A special case of such ruled surfaces is of importance and interest. If 
the ruled surface has the property of being tangent, along the rulings, toa 
moving plane which rolls around the surface, then it may be deformed by simple 
bending and flattened out into a plane. Such a process is called "development" 
of the surface. We call such surfaces 'developables" or "wrapped surfaces". 
(A sheet of paper can be wrapped around the two curves (u0) and (ul) to form 


the surface. These surfaces are also known as "convolutes".) 
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The tangent-plane condition can be established by showing that the tangent 
vector at a point on (u0), the tangent vector at a corresponding point on (ul), 
and the tangent vectors along the line joining these points, are all coplanar, We 
need to form the scalar triple product of these vectors, and show that it 
vanishes. ; 


The tangent vectors in question are (ud), (ul) > (ud) and (ul). 


First observe that, for a ruled surface, 

(uw), = (ul) - (u0). This tangent vector is independent of w, so 
that (ud) = (ul). Moreover, the vector is simply the line segment joining 
the two points, as might be expected. We can write, for the scalar triple 
product, 


(ud) 
(ul) = 0 
(ul) - (ud) 


where the notation represents the determinant of the matrix of the three (row) 
vectors. If the determinant vanishes for all values of u, the surface is 
developable. 


The preceding describes an analytical test to ensure that a ruled surface 
is developable. We shall now describe a construction that will enable us to 
define a ruled surface by means of two space curves. Suppose that the two 
Space curves are defined by vector functions of two different parameters, 

u and d 


$1 


ud 
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The scalar triple product is 


(u0) 


u 
@)g | = 0 
(P1) ~ (uo) 

If we consider u the independent variable, the equation enables us to find 
¢ for any value of u; this value of determines the point on (f1) which corre- 
sponds to a point on (u0), so that the line joining these points is coplanar with 
the tangent to (41) and the tangent to (u0). We allow u to vary, and obtain 
related $ values; these values of ¢ enable us to evaluate the components of the 
vector (#1). These components are the same as the components of the desired 


(ul) vector. 


Provided we remember that the (unspecified) functions of u and¢ are 
different, the symbolism 
(al) = Gl), ¢d = ¢ &) 
represents the statement that the [x y z] vector is the same for both. 
With the correspondence established between points on the two curves, 
we can write the equation for the developable surface, 
(uw) = (uljw + (uO) (1 - w). 
This is the ruled surface equation, but with a special relationship between 


curves (u0) and (ul). 
7.2 PLANE/SURFACE INTERSECTIONS 
The general surface equation can be cast in the form 


uw = uBw> where U and W are vector functions of u and w 
respectively, and where B is a square matrix describing the boundary curves. 


For example, we might be dealing with the first F-type surface equation, 


uw = - [-1 Fu Fu] 0 wd ul -1 
ow 00 Ol Fw 
lw 10 11 Flw 


in which these vectors and the matrix are explicit: 


74 SECTION VII 


Since uw is in reality a vector consisting of an x, a y and a z component, 


there are three B matrices which we can call BL By Bo. 


We wish to find the intersection of this surface with the plane 
ax + by + cz +d = 0. 
We can substitute x = UB Ww’, y= UBW" and z = UB Ww" into this equation, 
and write the result in the form 
U fs. + bBY + cB | we+d=o. 
It is permissible to interchange the order of multiplication from aU, bU, 
cU to Ua, Ub, Uc because a b and c are scalars. In this form, the sum 


(ap, + DB, + cB] = §, a square matrix function of u and w, and 


usw! = -d igs an equation in the two variables u and w. If w is 
assigned a fixed value, there results an equation in u which when solved will 
yield a point on the intersection curve of the surface with the plane, (provided 


of course such a point exists for the chosen value of w.) 


If the surface in question is a bi~cubic, the matrix S is no longer a 
function of the variables u and w, but consists of constant elements. In this 
case the above procedure reduces to the solution of a series of cubic equations 
in u, where the coefficients of the cubics are determined by successive fixed 


values of w. 


In any case, if the spacing of the w values is close, the old value of u 
just previously determined for a particular choice of w can appropriately be 
used as a first trial solution for the new value of w. Algorithms for the 
improvement of this initial trial value of u are not difficult to construct, and 


will not be discussed in detail. 


If the plane is given by, say, the equation 
x+d=0 
the solution procedure is unaltered. Not much simplification results. 
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SECTION VIII 
RATIONAL POLYNOMIAL FUNCTIONS 


8.1 BOUNDARY CURVES AND BLENDING FUNCTIONS 


Two kinds of curves have for many years been traditionally used in air- 
plane lines design — cubic polynomials, and conics. Unfortunately, each of 
these curve forms for itself has certain drawbacks. In the parametric form, 
for ordinary cubics, the entire shape of a curve segment is governed by end 
tangent vectors. Sometimes these end tangent vectors lead to unwanted hooks 
and bulges in the curve segments. On the other hand, conics, although more 
benignly behaved, cannot by their very nature yield curves with points of 
inflection. Yet such curves very often exist in aircraft shapes — as for in- 


stance in the case of wing fillets. 


Because of these short-comings, a new curve type has been developed. 
It is based upon rational polynomial functions. It contains both conics and 
ordinary cubics as special cases, and provides a great degree of generality 


and flexibility. 
We start by establishing the form of this function. 


Let v be a vector, so that for example v = [xy zl] orv = [xyl1] or 
v = [xl]. The first of these can be thought of as the vector (or matrix) of 
coordinates on a space curve; the second is the vector of coordinates for a 
plane curve, and the last is the vector of a single coordinate. Since this last 
vector yields the most general case, we shall begin with it, and show how one 
might evaluate a set of numbers in a matrix to define each of the parametric 


coordinates of a curve, 


The product of v and a variable scalar w is wv = [wx w]. Here both 


wx and w are cubic functions of a parameter, v, and obviously 
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This is the ratio of two cubic polynomials (hence the term "rational func- 
tion"). We can represent the two cubic polynomials by the matric equation 


we = [uwuzuly A. 
Since wy = [wx w], the matrix A must consist of four rows and two columns 
of constant coefficients. We now proceed to show how these numbers may be 
found so as to define a coordinate of a curve. , 
We shall be interested in the end-point coordinates of the curve atv = 0 
andv = 1, These coordinates are v. = [Xo 1] andv, = [x, 1] respect- 
ively. A tangent vector anywhere on the curve is clearly 


vii = [x' 0] 


where the prime mark means differentiation with respect to the parameter v. 
The tangent vectors atv = 0 andv = 1 are therefore y,. = [x,' 0] and 


v,' = [x," 0} respectively. 


Now differentiate both sides of 
wy = [u3u%ul} A and obtain 


(wv)' = (8u% 2u10) A. 


Substitution of u = 0 andu = 1 into these two expressions yields 


Wo Yo 0001 
wiv 11i1i1 
it . A. 
(Wy Vo)! 0010 


(w, v,)" 3210 
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The 4 x 4 matrix on the right has an inverse, and we may write 


“1 
0001 Wo % 
1111 wy, 
aes 11 

’ 

0010 (Wo Yo) 
t 

$210 (w, ¥;) 

221 1 Wo % 

= {78 3 -2 -1 WY, . 

t 1 

0 0 1 0 Wo Vo + Wo Yo 
' J 

1 0 0 0 wy, WtYy 


The square matrix inverse is constant and always the same, and reappears in 
the algebra so often that we shall henceforth call it the matrix M. 


The matrix equation can be factored and rewritten in the form 


A= M Wo 0 0 0 Fy 
0 Ww 0 0 % é 
: Wo 0 Wo 0 Yo 
0 w,' 0 Ww, y,' 


The right hand matrix of v's represents the desired end conditions on the 


curve. In our present case, it is of course a 4 x 2 matrix. 


The middle matrix is 4 x 4 and contains the four numbers [Wy ¥y Wo wi‘) . 
Any arbitrary set of four numbers inserted into this matrix will serve to define 
a unique pair of cubic functions of the parameter u, from which x can be found, 


by using, as we have said, the ratiox = = ° 


Instead of picking these four numbers arbitrarily, however, we shall 
impose further conditions on the curve until enough conditions are imposed 


1 ' 
to define [Wo W, Wo! W,') uniquely, 
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We begin by introducing desired second derivative vectors at the end- 


points; these vectors are clearly 


v," = [Xp" 0} and 


0 
Woe " 
eT {X, 0). 


(Incidentally, in the case of vectors v' = [x' y' 0], v'' = [x" y'' 0] if 
the determinant of the matrix e Al vanishes, the curve will have a 
x'' y" 
point of inflection at v. If the determinant is positive, the center of 
curvature will lie on the left as one proceeds along the curve; if the 
determinant is negative, the center of curvature lies on the right. If 
two curve segments have equal v' and v'' at a junction, they are contin- 
uous both in slope and curvature at such a junction. ) 


When we take second derivatives of both sides of 
wv = [u3u2ulj A, we obtain 
(wv)" = [6u 200] A, 
Atu = 0, this is 
t = 
(W¥)" = [0 2 0 0] M Wo Yo 
' 1 
Wo' Yo + Wo Vy 
t 3 
wi' Vv, + wv, 


" 
’ 


tos tt ! t " 
But (W) Vo) witty + 2wi vi + w Yo and solving for w, V) 


0 0 0 0 0 
Wo Vo" = (Wy Vy)" - Wo" fe 2 Wo Yo’ 
Now (Wo VQ)" = [-6 6 -4 -2] Wo Yo 
Wy My F 
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Furthermore, since in general 


(wv) = [(wx)” w'], 


the quantity w'' is the second component of the vector of (wv)"' and therefore is 


associated with the last column of the matrix [w.v 


0 and vy, are both 1, and the last components of Mie 


and Le are both zero. Hence 


But the last components of v 


w,"" = [-6 6 -4 =2] Wo 


We now can write, by combining results: 


yy tf . fie @ ud ~9 (ees Ao ok 7 2wiyt 
Wo Yo [-6 6 -4 -2] | ¥%9 % Vo) Lore 
®y Vv) oa) 
bg Ufa eS ; rd 
is (vy VQ) + Wy Vo 
PMS RE -s. A ' 
I, Vv, Vo) vant a 
= 6 wy (v mie - ws O ~ 2! iv, - Yo) -2w vy -2w 
Collecting, 
: Cee , Spas f z i ae ' 
Wo Yo Wo ( 2 De AO Ms Vo) vied 


wb Foo cal rh Elo ty oo 
+ (-2 V9") + wy (-2 (Vv, V))- 
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We now restore this last expression to matrix form: 
W t ' t 
Wo%o = [W WW Wy ] |-4 Y% 


ag ‘aa t 
.6 (v, VQ) 2v, 


= ' 
2 Y 
~2 (v, - ¥y) 


In our present case, the matrix on the right consists of a column of numbers 
and a column of zeros. Hence the column of zeros can be discarded, and the 


result is a 4.x 1 matrix. On the left, Wo var ig a scalar. 


Similarly, we can find by analogous algebraic procedures that 


te t ' = ' 
wi, = [Wp wy Wo Ww, ] 6 Vo v,) +2 Vo 
4 v;' 
2 Vy - v,) 
= ' 
2v, 


Then, writing a matrix equation, we have, so far, 


[Wy Vo" W,¥,"1 = [WW Wo'W,') (P |Q] 


where [P|Q] represents a 4 x 2 matrix consisting of the separate 4 x 1 matrices 


for Wo Yo" and Wy, vy" written side by side as columns. 


We now introduce another condition. Let it be required that the curve 
pass through the point is [x, 1] whenu = + - (This value of u is of 
course arbitrary. ) 


This condition leads to 


i 
1s 8 (1248) M w, Vv, 
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By algebraic manipulations similar to the preceding, we can rearrange the 


equation to read 
= t ! ' 
as [Wy ¥, Wy w,'] Avot 
a ’ 
4 vy % 
i) 
mel 


This is an equation in which Ye hag two components, x, andl. The right hand 
matrix is a4x2. Call it the R matrix. Then we can adjoin these matrices, 
to obtain 
t 
[Wo Yo i wiv 


" 8vi) = [WW Wo Wy"! [PQ R]. 


1 
Now [PQ R] represents a 4 x 4 matrix; P and Q are each 4x1 matrices, but 
Risa4 x 2 matrix. We next transfer Wo VQ" and Ww v," to the right hand 
side, obtaining 


= ' ' - 
[09 9 8v] [Wy Ww, Wo wy ' [PQR] Vo" 0 0 0 


The right hand matrix is now a 4x4. Provided its determinant does not vanish, 


it has an inverse, and 


[> Wy Wo w,'] = [0]0 sv] Ss 


where S is the 4 x 4 inverse of the matrix. 


Now that [ Wy 3 
pletely defined, since the rational function 


Wo and w,' ] have been evaluated, the curve is com- 


is completely defined. 
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8.2 PLANE CURVES 


With some loss of generality and flexibility, we can have v = [x y 1], 
a plane curve. We shall next show that for an appropriate choice of 


[Wy Wy, Wo! and W,'] , the curve reduces to a conic. 


We have the equation 
wy = [wx wy w] 
= [ w uu ] A. 
In this case, A is a 4x 3 matrix. Now if the top row of this matrix is 
[0 0 0], the equation reduces to 
wy = [ wv uljA 
when the top row of A has been omitted. A is now a 3 x 3 matrix, and it is 


possible to show that this equation is a parametric form for the general conic, 
expressed as a quadratic rational function. 


For the top row of A, we have the vector equation 


[2-211] [wv = [00 0] 


Expanding: 


f t t t a 
2 Wy Vo 2w, vy) +Wo Yo t Wo %y + Wy Y, +, Vy [0 0 0} 


Collecting: 


' ex, t ' ' a 
Wo (29 + V9') + W, (-2¥) + ¥)') + Wo (Vo) + (vy) [0 0 0] 


In matrix form, 
rn ’ t ’ 
[0 0 07 [Wy Ww, wi") 2v)+Vo 
ei t 
2v,+V) 
%9 
Vv. 


1 
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The matrix isa4dx 3. 


Now we can adjoin a column to the matrix and an element to the vector, 


and write 


= y * ae a) 
fo 00 wo = [Wy Wy Wo wii) 


then 


If the indicated inverse exists, then a solution can be obtained in terms of w 


Furthermore, Wo can be set equal to 1 arbitrarily. 
The matrix has an inverse in case the determinant 


-2Vv pot , 
2v, +, 4% 0 


As atest, construct a conic with end conditions 


, r001] 
% | 1 | 
| 
1 


| 1 
vy Pia] 
bg 100 
' 
4 Vv. o————_———__> \ 
v,! } 100! 0 0 
1 u | 


vy | 0 


al 


0° 
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2v,t+v'! [002] 4,100] 
“2v, tv)! [-2-2-2] +[10 0] 
Vo 001 
Vv 111 
14 L 


We first test to see whether the determinant vanishes: 
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0 2 
-2 -2 
Oo il 
pees E 


-1 -2 =2 
~1 -2 
0 oO 1 = - = -l. Hence the augmented matrix 
at a 
1 i 1 
will have an inverse. 
The matrixis | 1 0 2 1 and its inverseis | 0 1 0 2] 
-1-2-2 0 0-1 =1 =1 
001 0 001 0 
111 0) 1-1-2 -23 
then 
+ t See 
[1 w, wo’ wt] = [0 0 0 1) 0102 
0 -1 -1 -1 
0010 
el 
h ' ! = eh = = 
whence [1 Ww, Wo wid [1 -1 -2 -2] 
The conic equation is 
wy = fusu2ulj | 2-2 11 1000 meee 
-3 3-2-1 0-1 0 0 111 
0 01 0 -2 01 0 1 0 0 
100 0 0-2 ~-1 10 0 
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= {wu2ul] fo 0 1-1 001 
1-1-2 1 111 
2010 10 0 
1000 100 

= fueu2u1y [0 0 0] = futnay [-2-1 0 
-2-1 0 1 0-2 
1 0-2 001 
001 


[wx wy w] = [-2u? + u | -u? ] -2u +1] 


x =u 
au =x" 
YS rae vet ee Ft The curve is hyperbolic, with an 
asymptote at1- 2x = 0, x = +. 


It is always possible in all of the foregoing to set Wo = 1. This is 
because all equations are homogeneous. It is never possible for Wo = 0, 
since this leads to certain degenerate cases. 

We remark in passing that when 

[Wo Wy Wo wi] = [11 0 0] 
the equation reduces to the ordinary parametric cubic, given by 


v= [we u2uiy M fv 


and w is constant and equal to 1. Hence the rational polynomial functions 
contain as special cases all conics, ordinary cubics, and of course therefore 
straight lines and circles. 
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Their use as boundary curves for surface patches is obvious. They 
maintain tangent vector continuity between adjacent patches; indeed, if the 
F 0 and Fy functions are constructed as rational functions, we can establish 
the F function 


This function has the end conditions 


Yo 001 
vy 111 
Yo = 10 0 ; 
‘i 100 
00 0 
vs 00 0 
Since v)" = Avy! = Of1 0 0] andv," = dv! = 0f1 0 0}, the 


curve has a point of inflection at u= 0 andu=1. Hence its use insures 
curvature continuity across boundaries between patches, provided of course 


the boundary curves have similar curvature continuity at patch corners. 
The cure is symmetric, Furthermore, we can put 
Fo(u) = 1 - F,(u) 


and obtain directly the Fy function, another cubic rational function, with 
similar properties to F,- 


8.3 AN EXAMPLE 


We shall work out the equation for the F, (u) blending function with the 
customary stipulations that 


F,@) = 1, Fi (0) = F,'(0) = F,'(Q) = 0, 


and with the two additional stipulations that F,"(0) = FQ) = 0 as well, 
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This blending function will give both slope and curvature continuity across the 
common boundary between two contiguous patches, The end-conditions are, 


forv = [ F(a) L]y 


7 [f- =) 
[ %o le Z 
i 
| 
y t 1 “| 
vt 10 O: 
pi | 9 0 | 
py I 
Vou 0 0 
L“i J Le 
The matrix [P Q R] - aes 0 0 0 | 
0 pt 0 
vy 0 
jo 90 0 0 
0 o 0 0 


obtained by direct substitution in the given form. 


Its inverse is 0 -1 oO -4 = § 
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1| for symmetry. Then, since u = -=~, 


t Oe ae 
[Wy Ww w ws [0 8 ves 


[0 0 4 8] 0-1 O = 


a 
Oo 
i 
oD oO fF 


= [8 8 -24 24]. 


Now if Wy = 1, instead of 8, the equation becomes 


[Wo W, Wo wi") = [11 -3 3] 


we have 


and substituting the values of [Wo Wy Wo w,' I 


ye ee eas a Ct 2 ja % 
-3 3 -2 <1 1 1 0 3 
0 oO 41 0 0 -3 0 =3 
1 0 oO 0 3. 3 0 dil 
Finally, wv = [wx w] = [ue uw ul] 1 Of 
0 3 
0 -3 
0 1 


[wx w] = fue | 3ue - 3u 41] 
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3 
Hence x = —e- = ——t—_——_ = F,(u) as required. The other 
w 1 
suv = Sut] 
F, function is 
: 3 
Fy(u) = 1-F,(u) = 1- = 
su -3u + 1 


ee oe ee 


gue - su 41 


This is seen to be a rational cubic function also. 
8.4 PLANE CURVE THROUGH A POINT 


The plane curve vector isv = [x y 1], wv = [wx wy w]. Here 
the polynomial denominators inx = = andy = = are both the same. 


As before, the end conditions on the curve are contained in the 


matrix Vo a 


We wish to cause the curve to pass through some arbitrary point Yo 
(commonly called a "shoulder point") and it will turn out that we shall also 
be free to choose some arbitrary slope at this point. It is important to 
distinguish between the term "slope" and "tangent vector". The slope of a 
curve implies that the direction of the tangent vector is known, but the magnitud 


of the vector is not under our control. 


We begin by assuming some value of the parameter u to correspond to 


Vor For purposes of illustration, let u = + at this point. 
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Then, from wv = tua’ ui] A, we write 


1 

Yor fl 24 81M Wo Yo 
bes eh 

' ' 

Wo' %y + Wy % 


t ’ 
Ww Y% + Wi vs 


Observe that we have arbitrarily set vo 1. This is harmless, since 
the equation is, as we have observed, homogeneous. 


By multiplying the matrices, combining, and collecting terms as we 
have done before, we achieve the result 
= ! 1 t 
Sy. [W, Ww, ¥% wid Avy +V, 


S ' 
4v, % 


The matrix on the right is a 4 x 3 matrix; to make it square, so that it 
can have an inverse, we need an additional column. This column can be pro- 
vided by a scaler equation, and the slope relationship will furnish this equa- 


tion. 

We first find an expression for the tangent vector at v et Differentiating, 
we obtain as usual, 
uz 2u10] A 


(wy)! = 
w' = [ $u” 2u1 0] M Wo 
' Ww 

t 

Wo 

w.' 
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This last equation comes from the equation for (wv)' by an argument that 


we have used before; that since wv = [wx wy w], (wv)'’= [ (wx) (wy) wi]. 
Hence w! corresponds to the last column of [w, vi 1 
bee 
ie ' 
(Wy Vo) 
' 
Oy v,) 


which is [| Wo 


1 : age Se 
Now at a we have already set u = > «We make this substitution, and 


obtain, from 


wv! = (wv)! - wiv 

ae 6 6 -1 1 [w, (v_-v_) 

ee te 0*'o Ve 
Ww, Vv, -v) 


|, aa: tt = P ' 
| %o Vo ¥) + NG 5 


t - A t 
fl (v, %) + Ww Vy 


When we perform the indicated multiplications, and then collect results 


and restore to matrix form, we have 
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We are now ready to introduce the slope condition. We could write 


x . x , but this would yield awkward results when the slope became very 

great and approached an infinite slope. Instead, we choose two numbers a 

and b so that ax’ = by'. It is obvious that these numbers can very appropri- 

ately bea = sin @ , b = cos @ where @ ig the slope angle. Then, for Vo" 
' 


O= by, ~ ax 


This is a scalar equation. 


Now xO is the first component of vo and y, is the second component 
of vo Hence these quantities correspond to the first and second columns of 
the matrix in the equation for ve respectively. We can write this out in 
detail: 

a ' ' os ay ft 

0 [Wy Wy Wo Wy, ] b 6%, ¥) Yo 
‘ie = t 
6, -¥)-¥) 


FeV 
Ig 7 
bat ~6(x)-x,) - x, = [> WwW, Wo w,' ] p 
6(x,-x,) - x,' 
r 
x, 7 Xp 8 
x7 % 


where the brackets on the right enclose the resulting 4 x 1 matrix (or column 


vector), 


The factor 4 (of 4 v,') obviously drops out of the equation. 
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We now adjoin this to the equation for YS Again we can drop the 8 (from 
8 Vis and obtain 


[Vv 0] = [Wo Wy Wo w,'] 4vy+vo' p 
4v,-% q 

Yo r 
“vy 8 


The vector on the left consists of four components: (x, 


matrix on the right is a 4 x 4 matrix. If it has an inverse, S, then we can 


a 1 0]; the 


t t 
solve for [Wy Ww) ¥y ] by the equation 


1 ' ec 
[Wo Ww) Wo Ww, ] EY, 0) S. 


These values of the w vector cause the curve to satisfy the desired 


conditions, 
8.5 SECOND DERIVATIVE VECTORS 
We have already discussed rational functions for 
v= (xX 1] (and ofcourse [y 1] and [z 1].) 
In particular, we showed that curves based upon these functions can 
usually have arbitrary first and second derivative vectors at the end-points, 


and in addition can be caused to pass through some shoulder point Yoo also 
arbitrarily chosen. 


When the vector v = [x y 1], the complete generality of the resulting 
curves is somewhat curtailed. We shall investigate the conditions under which 


such a plane curve can satisfy end conditions including second derivatives. 
We have already obtained an expression for the vector 
" w = t t 
[Wy %o" W, ¥,") = [% Wy, Wo! W,') {P| Q}. 


Before, the vectors Wo Vo" and Wy v," were actually scalars, since they came 


from v" = {x" 0]. But since v' = [x" y" 0], they are each 2~component 
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vectors, and their combination makes a 4-component vector. Similarly P and 
Q are now each a 4 x 2 matrix, and their combination is a 4 x 4 matrix. 


We carry w, v,." and Wy v." across the equal sign, and obtain on the left 


0-0 1 
the null vector: 
' - 

[0 0 0 0) = [w, W, % Ww") [Pa] Vo" | 0 0 

wt 

0 0 Y, 

0.0 0 0 

0 0 0 0 


The matrix on the right is 4x4. Now the condition that must hold, in 
order for there to be a solution for [> Wy Wo w,' ], is that this matrix 
must be singular; the determinant of this matrix must vanish. 

This last remark tells us that % 
arbitrarily. However, it is always possible to make the determinant of the 
matrix vanish by the adjustment of any one of the four components of Vo" and 
v,''. Thus if one of the four components is the number a, we can expand the 


1 
determinant in such a way as to obtain the equation. 


" and v," cannot be chosen entirely 


k, a+ k, = 0, from which a can be found. 
Suppose the matrix is, or has been caused to be, singular. Then, if 


[0 0 0 0} = [Ww w w,' wi'] §S, 


1 °0 1 ] 
we make it non-singular by an appropriate modification. In some cases, this 
might consist in adding 1 to an element in the top row of S. As an illustration, 
we might have 


2 t t 
[9 0 0 wo] [Wp W, Wo wi’) S + 


coco fo 6 
o oc 8 S&S 
— ~~ 
le ~ er 
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It can be seen that this modification is still a valid equation; if the 
modified S matrix now has an inverse, we can immediately obtain it and solve 
for the w vector. 


As an illustration, consider the end conditions 


% 001 
vy 111 
Vo! 100 
y,' : 1 0 0 
VQ" 100 
vy, a 0 0 


We plan to adjust a in v," until the matrix is singular, 


We require first the matrix [ P | Q], which is given by 


-4 v,' 6 (% - v,) +2 Vo! 

6 (v, - Vy) -2 v, 4 v, 

-2v,! 2 (Vy -¥,) ; 
~2 (v, - V9) -2 vy! 


Substitution of the end conditions gives the 4 x 4 matrix: 


-4 0 [-6 -6] + [2 0] -4 0 -4 -6 
[6 6) -[2 0] 4 0 4 6 4 0 
2 0 -2 0 =2 2 O -2 2 
“2 -2 2 0 2 -2 -2 0 


In passing, we note that if our end conditions had been Vo" = [0 0 0} 
and if also v," = [0 0 0}, the resulting matrix would be singular, because 
the first and third columns of { P| Q] are identical, However, this is not 


our present case. 
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We now subtract the matrix 


we 
Vo 0 0 1 0 
Ad 
0 0 vy 7 0 0 | 
00;0 0 0 0 
0 0;0 0 0 0 


-5 0 ~4 -6 

4 6  (4-a) 0 
“2 0 -2 -2 
-2 -2 “2 0 
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By a series of reductions accomplished by multiplying rows of the matrix 


and additions (or subtractions) of rows to remove elements, we can obtain the 


deter minant 


(a-1)-1]/ =0 


1 -1 


This implies 


a= 2, 


This is the value of a that makes the matrix singular. 


The singular matrix is 


-2 


We make it non-singular by adding 1 to the top left element. 


is now 
4 0 -4 -6 
4 6 2 0 
-2 0 -2 =2 
—2 -2 -2 0 


-6 
0 
~2 
0 


The matrix 
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and corresponds to the vector [ Wo 0 0 0] onthe left of the equation. 


The inverse of this matrix is 


to: . | 
3 2 1 -6 3 
-l 0 3 -1} =R 
-1 -1 3. -3 | 
-1 0 2 0 
Finally, [Wp Wy Wo! w,'] = [Wo 0 0 0} R. 


If we arbitrarily set w. = 2, then the required solution is just the top 


0 
row of R, or 


= [21 6 3}. 


With these numbers known, the curve equation is completely defined. We 


obtain it by substition in the canonical form: 


WV = [wu ul] M [wv 


[0 0 -6] + [2 0 0] 


= [uuwul) le = ee ee See ee 
23. 8 48-21 a 
0 0 1 Oj {2 0 -6 
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2 
pueue ul] 


2 
0 


mH 


coco 2 0 


~1 


-6 


iw) 
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The separate equations for the x and y coordinates of the curve can be 


exhibited: 


aye + 6a 


9 
Suh 4 6u" = 6u + 2 


3 
u 


3 2 
4u  -5u +2u 


~ 6u +2 


In the foregoing, certain matrices have occurred. These matrices are 


significant ones, and can be written as transformations of the common 


a. 1 
matrix [ v) ] 


as follows: 


For the conic condition matrix, 


0 


-2 
mf 


r2 45] 
Vv 2NG 


+v_' 
a4 
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For theu = > 


2 


For the P and Q matrices associated with w 


7 1 
4 
~ Fo tees, 1 
6Vv 6 Vv, zy, 
-? ' 
%9 
~2 
eles ass 
6 6 2M 
io vi +ey 
4v_' 
2v,-2V¥ 
-2 v,' 


0) 


4 0 
0 4 
1 0 
0 = -1 


0 -4 
6 0 
0 -2 
2 0 
6 2 
0 0 
-2 0 
0 0 


Me shoulder point condition, 


4 


-2 


tw 
0 %o 


and w_ Was 
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